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ANALYSIS OF THE MARKET WEIGHTS UNDER THE 
VOLATILITY-STABILIZED MARKET MODELS 

By Soumik Pal* 

University of Washington 

We derive the joint density of market weights, at fixed times 
and suitable stopping times, of the Volatility-stabilized market mod- 
els introduced by Fernholz & Karatzas in 2005. The argument rests 
on computing the exit density of a collection of independent Bessel- 
square processes of possibly different dimensions from the unit sim- 
plex. We show that the law of the market weights is the same as 
that of the multi-allele Wright-Fisher diffusion model, well-known in 
population genetics. Thus, as a side result we furnish a novel proof 
of the transition density function of the Wright-Fisher model which 
was originally derived by Griffiths by bi-orthogonal series expansion. 



1. Introduction. The family of multidimensional diffusions models named 
Volatility-stabilized market (VSM) was introduced by Fernholz and Karatzas 
[13] as a toy- model that nevertheless reflects some of the traits of a real- world 
equity market. We refer the reader to an excellent survey article by the same 
authors [14]. These models reflect the fact that in real markets the smaller 
stocks tend to have a greater volatility and a greater rate of growth than 
the larger ones. 

The mathematical description of the model involves a vector-valued con- 
tinuous stochastic process X{t) = {Xi{t), X2{t), . . . ,Xn{t)) where every co- 
ordinate takes nonnegative values. Their dynamics are determined by the 
following stochastic differential equation with a single nonnegative parame- 
ter 6. For i = 1, 2, . . . , n, we have 



(1) 



dX.it) = -S{t)dt + ^/X,it)Sit)dWi{t), Sit)=X,it) + ...+Xn{t). 

The initial vector, X{0), is a point in the positive quadrant of M."", which we 
will denote by M""*". Here {Wi, W2, ■ ■ ■ , Wn) is an n-dimensional Brownian 
motion. The original article by Fernholz and Karatzas [13] parametrizes the 
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model hy a = 5 — 1 which is assumed to be nonnegative. Our analysis will 
consider a more general class of models where the scalar 6 is replaced by 
a vector ((5i, . . . , (5„) of nonnegative coordinates with only the restriction 

Er=i s^ > 1- 

The intuition behind such a modeling becomes clear from the following 
consideration. Define the vector of market weights 

(2) i^i= ' i = l,2,...,n. 

From an economic viewpoint, the market weights are a measure of the in- 
fluence that the ith company exerts on the entire market and have been 
studied extensively in the literature. For example see articles by Hashemi 
[19], Ijiri and Simon [20], Jovanovic [21], and Simon and Bonini [27]. For 
a probabilistic study in the context of another interacting market model, 
please see the article by Chatterjee & Pal [6]. 

One can alternatively express SDE (1) (see [13]) by writing 

dlogXi{t) = P^dt + -^l=dWi{t), i = 1, 2, . . . , n, 

which makes some of the features of the model immediate and visually ap- 
pealing. The smaller /Xj is, the greater is the drift and the fluctuation of 
log Xi . This is the primary empirical observation that the model is designed 
to capture. 

In this article we answer one of the questions left open in the articles 
[13] and [14]: how to describe the behavior of the vector of random market 
weights (/ii, . . . ,Hn){t) under the law of the VSM model? Similar problems 
have been studied by Irina Goia in her thesis [16]. Please see her thesis 
for a discussion on the relationship of these models with the CIR models 
in mathematical finance and their relevance in the bigger picture of the 
stochastic portfolio theory. 

As a natural culmination of the theory we develop in this article we con- 
sider the following generalization of the VSM models. 

Definition 1. For any n nonnegative parameters (Ji, . . . , consider 
the solution of the stochastic differential equation 

(3) dlogX^{t) = ^^dt + ^l=dWi{t), i = l,2,...,n. 

We call the unique in law solution of the above equation the VSM model 
with parameters . . . , 6n) and denote it by V((5i, . . . , (5„). 
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As mentioned in [13], the uniqueness in law of the above SDE is guaran- 
teed by results in the theory of degenerate stochastic differential equations 
as developed by Bass &: Perkins in [2]. 

A crucial observation made in [13] in analyzing the VSM model is their 
connection with the Bessel-square (BESQ) processes. Given a solution of 
SDE (1), one can construct n independent BESQ processes of dimension 26, 
(say) Zi, Z2, ■ ■ ■ , Zn, such that the solution X is linked with Z = {Z\, . . . , Zn) 
by an appropriate time-change. Explicitly 

Xi{t) = Zi {A{t)) , 0<t<oo, i = l,2,...,n, 

(4) 1 /■* 

A{t) = - S{u)du, S{u) = Xi{u) + ... + Xn{u). 
4 Jo 

A straightforward generalization of the analysis of Fernholz & Karatzas 
shows that a weak solution of the system in (3) can be obtained by the 
following mechanism. Given a solution X of V{5i, . . . ,6n), there exist pro- 
cesses Zi,...,Zn, which are independent BESQ processes of dimensions 
(2(5i, . . . , 2(5„), such that the time-change relation described in (4) contin- 
ues to hold. 

We have the following results. 

Proposition 1. Let X = {Xi, . . . have the law V{6i, . . . ,dn) as in 
(3) with initial Xj(0) = Xi > for every i. Suppose 

n 

6i > 0, for all i, and d = 6i > 1. 

i=l 

Let S{t) denote the total sum process Xi{t) + X2{t) + . . . + Xn{t). Let ?a be 
the stopping time 

n 

c,a = inf > : S{t) = a} , s := Xi < a. 

1=1 

Then, the joint density of the market weights fi = {fii, . . . , jin) at the stopping 
time c^a is given by the following expression: 

^x[y) = 1 - s a) 2^ — — — — - (1 + s a) 
^-^ m!l (m + a) 

m=0 ^ ' 

(5) X (fci.'^.fejn^^^/'')'' ^''■(^'^+^)' 

fe>0: fei-|-...-|-fen=ni ^ «=1 

n 

Vi ^ 0, for all i and yi = 1. 

i=l 
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Here, k + 5 denotes the vector (fci + 5i, . . . , /c„ + (5„), and Dir{y;^) is the 
density of the Dirichlet distribution with parameter 7 given by 

T-rn p/ \ _n n 
(6) I)^,(y;^) = il^^l^ JJy7«-\ y,>0, ^yi = l. 



As mentioned in the abstract, the analysis requires us to compute the 
exit density of a cohection of independent BESQ processes of dimensions 
. . . , 5n) which might be of independent interest. 

Proposition 2. Suppose n > 3 and let Z = (Zi, . . . , Zn) be independent 
BESQ processes of dimensions {61, ... , On) where 



> for all i, and Oq = 6i > 2. 



i=l 



Assume that initially Z{0) = z = (zi, . . . , z^), where each Zi'>Q and ■- 
Y17=i < 1- Consider the stopping time ai given by 

f7i=inf{t:C(t)> 1}, at) = Zi{t) + ... + Zn{t). 

Then, the density of {Zi, Z2, ■ ■ ■ , Z„)(cji) is given by 

/ N /I ^ r(2m + 6*0/2) c- ^-2m.-eo/2 

^-^ m!i m + On 1) 

(7) X Yi L."'.feJn^^^^'^(^'^+^/2)' 

fc>0: fci+...+fc„=m ^ 1=1 

n 

Ui ^ 0, for all i and yi = 1. 



Here, k + 9/2 denotes the vector (ki + 9i/2, . . . ,kn + 9n/2). 

Since each 9i is assumed to be strictly positive, the above expression is 
also the exit density of the Z process from the unit simplex {x G M" : Xi > 

0, Er=i^i<i}- 



A deeper analysis can be done by noting, as we will show in Section 3 
that the distribution of market weights under the VSM model is nothing but 
the multi-allele Wright-Fisher diffusion model studied in population genet- 
ics. A short introduction to this well-known and important model has been 
provided in Subsection 1.2. 
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Proposition 3. The process of market weights {fii , . . . , fin) under V{5i , . . 
is itself a diffusion independent of the total sum process S. Its law is the 
same as that of a multi-allele Wright-Fisher model with mutation parame- 
ters . . . 

Under the additional assumption that each 5i is strictly positive, the unique 
reversible invariant probability law for the market weights underV{6i, . . . , 6n) 
is given by the multivariate Dirichlet distribution with parameters . . . , 

Finally, we prove a transition density formula for the market weights. 
Since we show that the market weights have the same law as the Wright- 
Fisher diffusions, it follows that this is the same as the transition density for 
the Wright-Fisher model which was originally derived by Griffiths in '79 [18]. 
Also see Griffiths [17]. Our proof is novel and follows easily from Proposition 
1 and suitably changing time. 

Proposition 4. Let p{t, y) denote the transition density from an ini- 
tial point ^ to a final point y of the market weights under the VSM model 
which satisfies the same assumptions as in Proposition 1. Then p{t , ^ , y) is 
given by the formula 

^-^ mil im + a) 
(8) X L.^'.A: )nte^''^'^(y;^ + '^)' 

fc>0:fci + ...+fc„=m ^ i=\ 

n n 

(i>0, >0, = ^y, = 1. 

i=l i=l 

The coefficients brn{-) can be expressed by the Laplace transform formula 
which holds for all positive p: 

^p-lg-(m+7)p _ g-p) + g-p) -2m-d ^ ^ = 0, 1, 2, . . . . 

Here 7 = {d-l)/2. 



Remarks. 

(i) Tavare [28] has a different proof of the above formula for the Wright- 
Fisher model where the coefficients bm{t) are themselves linked to transition 
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probabilities of a pure death process in Z"*" U {00} . Our formula above estab- 
lishes a Laplace transform representation of the same probabilities, which 
might have an interest. 

(ii) The transition density function for the Wright-Fisher model, as de- 
rived by Griffiths, has exactly the same form for all nonnegative values of 
{5i, . . . ,6n)- It should be possible by extending our methods to get rid of 
assumptions on the parameters. But it is not immediate and requires some 
more work. We do not pursue it here since the VSM models naturally as- 
sume that XlILi > which corresponds to the fact that the entire equity 
market never hits zero. 

(iii) There is an interest in determining whether the market weights in equi- 
librium exhibit power-law decay. That is to say, the ith largest market weight 
fii is proportional to , for some positive 7. This is empirically observed, 
and can be proved in the case of certain models. See Chatterjee & Pal [6] 
for more motivation and references and some results involving the Poisson- 
Dirichlet families of point processes with parameters {a, 0) where this in- 
deed takes place. However, there does not appear to be such a possibility for 
the VSM models. The finite-dimensional invariant distribution have been 
identified in Proposition 3 to be the Dirichlet distributions. Under standard 
Poisson convergence assumptions, the point process of the order statistics of 
Dirichlet distributions converge to the Poisson-Dirichlet's with parameters 
(0, for some positive /?, which do not exhibit power-law decay. 

The paper is arranged as follows. The next subsection describes the multi- 
allele Wright-Fisher models and their limiting measure- valued diffusion, the 
Fleming- Viot model. In Section 2, we provide proofs of Proposition 1 and 2. 
This is achieved by defining a multidimensional functional transformation, 
akin to the Kelvin transform for the Laplacian, that utilizes inversion with 
respect to the unit simplex. In Section 3 we establish the fact that the 
process of market weights under the VSM model is actually the Wright- 
Fisher model. The analysis is slightly generalized to include the Fleming- 
Viot models which shows the large n behavior of the market weights. In 
Section 3.1, we consider the practical situation when one considers not the 
entire vector of market weights, but only a subset of it. It is shown that 
there is a recursive property of the VSM models which allows us to handle 
it effectively. Finally in Section 4 we establish Proposition 4 as a corollary 
of the previous results. 

1.1. A collection of notations.. This article sometimes requires notations 
that refer to similar, and yet different, objects. To help the reader avoid 
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confusion, we list most of the notations used repeatedly in the following 
sections. 

The VSM processes will be denoted throughout hy X = {Xi, . . . , Xn) 
while the BESQ processes will be Z = (Zi, . . . , Their dimensions will 
be the vectors 6 and 6 respectively. The sum processes will be 5 = Xi 
and C = Y17=i with corresponding dimensions 

n n 

(10) d = ^5i and 6^ = ^61. 

i=l i=l 

The stopping times ?a and aa denote the random hitting time of level a 
by the processes S and ^ respectively. Sometimes it will be convenient to 
consider the following transformation of the parameter 0: 

n 

(11) u,:=9j2-l, uo = ^Ui. 

i=l 

The closed positive quadrant in n-dimensions will be denoted by M""*". We 
denote the n-dimensional closed unit simplex by 

(12) 

S = < X = (xi, . . . , Xn) ■ > 0, for alH = 1, . . . , n, and Xi < 1> . 



i=l 



The oblique boundary of the unit simplex will be denoted by 
(13) 

6 = < X = (xi, . . . , Xn) ■ Xj > 0, for alH = 1, . . . , n, and Xj = 1 > . 

For any two vectors a = (ai, . . . , an) and b = (61, ... , bn) we will use the 
following notation: 



(14) Sa = Y,ai, a' = ll4\ a! = 11' 



a" = lla^', a\ = J_J_a,:!. 

i=l i=l i=l 



1.2. A brief description of various models. In this subsection we describe 
the various stochastic processes which are all linked to the VSM models. 

(i) The Bessel-square processes. A comprehensive treatment of BESQ 
processes can be found in the book by Revuz & Yor [25]. This family of 



imsart-aap ver. 2010/08/03 file: vsm_final.tex date: August 5, 2010 



8 



one dimensional diffusions is indexed by a single nonnegative real parame- 
ter 6 (called the dimension) and are solutions of the stochastic differential 
equations 

(15) Z{t) = x + 2[ y/\Z{s)\d(3{s) +et, x>0, t>0, 

Jo 

where /? is a one dimensional standard Brownian motion. We denote the law 
of this process by BESQ^.. It can be shown that the above SDE admits a 
unique strong solution which remains nonnegative throughout time. 

For = 1, 2, 3, 4, . . ., however the same process law can be obtained from 
another perspective. It is well-known that in dimension 9 = 1,2,3,4,..., 
the BESQ process has the same law as that of the square of the Euclidean 
norm of Brownian motion in dimension 6. The case = is unique. BESQ 
process for dimension zero is a nonnegative martingale which is a diffusion 
approximation to the process of the size of the surviving population of a 
critical Galton- Watson Branching process. 

The applications of BESQ processes, and especially of the derived Bessel 
processes, are too numerous to recount. To point toward the tip of the ice- 
berg, we mention such diverse areas as: (i) Branching process theory and 
superprocesses. See Etheridge [10]. (ii) Brownian path decomposition and 
excursion theory. See the book by Revuz & Yor [25, chapter XII]. (iii) 
Levy processes, see the article by Carmona- Petit- Yor [5]. (iv) Local times 
of Markov processes and Dynkin's isomorphism. See Eisenbaum [9], Pitman 
[23], and Werner [30]. (v) Mathematical finance. See Cox- Ingersoll- Ross [7], 
Geman & Yor [15]. (vi) Random matrices: see Bru [4] and Konig & O'Connell 
[22]. 

(ii) The Wright-Fisher diffusions. The Wright-Fisher diffusion model 
(see, e.g., Ethier & Kurtz [11, page 432]) arises as the diffusion approx- 
imation of the Wright-Fisher Markov chain model as the population size 
goes to infinity. A good source for an introduction to the Biology and the 
mathematics of these models is Chapter 1 in the book by Durrett [8] . 

For the purpose of this article it is a family of diffusions with state space 
& and parametrized by a vector (5i, . . . ,6n) of nonnegative entries. They 
are the solutions of the stochastic differential equations: 

1 

(16) dJ{t) = -{5i-dJ{t))dt + a{J)dl3{t), d = ^5^. 

1=1 

Here /3 is a standard multidimensional Brownian motion and the diffusion 
matrix a is given by 

(17) = {Hi = j} - ^XiXj) , l<ij<n. 
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The law of this process wih be denoted by J{6i, . . . , (5„). 

Sometimes in the hterature this process is identified by its Markov gen- 
erator: 



For the case of ?7, = 2, the first coordinate of the Wright-Fisher diffusion 
is also known as the Jacobi diffusion; see the article by Warren & Yor [29]. 
Hence, sometimes the general class is also referred to as the multidimensional 
Jacobi diffusions. For example, see Goia [16]. 

It is known that, for any n € N, and any strictly positive 5i,...,6n, 
the Dirichlet distribution Dir((5i, . . . , (5„) is the unique reversible invariant 
measure for the Wright-Fisher model J{Si, . . . , (5„). See Lemma 4.1 of [11]. 

(iii)The Fleming- Viot diffusions. The large n limit of the Wright-Fisher 
diffusions is the family of measure-valued diffusions that go by the name 
of Fleming- Viot processes. See the survey by Ethier & Kurtz [12]. These 
diffusions take values in the set of all probability measures on an underly- 
ing space and can be parametrized by a linear operator. The Fleming- Viot 
process and the Dawson- Watanabe superprocesses are probably the most 
important families of measure-valued diffusions studied in probability. For 
an introduction to the rich literature in this area, please see the book by 
Etheridge [10]. 

We hardly need the general theory. In fact, the family of Fleming- Viot pro- 
cess we use in this article has no spatial component. Let B be any Lebesgue 
measurable subset of [0, oo) whose Lebesgue measure is for some ^ 0. 
By a Fleming- Viot process we will refer to a stochastic process which at 
any point of time takes value in the metric space of V{B), the set of all 
probability measures supported on B under the Prokhorov metric of weak 
convergence. This process, say v, is defined by the following property: for 
any n = 1, 2, . . ., and any partition of B into disjoint Lebesgue measurable 
sets Ai, . . . , An of Lebesgue measure 6i, . . . ,6n respectively, where 9i > 
and ^0 = S"=i ^i; the law of the derived process 



is distributed as J{9i/2, . . . ,6n/2). We will denote the law of the process 
by FY{B). We will construct such a process later in the text which will 
prove its existence. That it is uniquely defined by the above specification is 
clear. 
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2. Description of the market weights under the VSM models. 

Consider n nonnegative parameters (^i, . . . ,6n), and n independent BESQ 
processes {Zi, . . . , Z„) where the dimension of Zi is 6i and the assumptions of 
Proposition 2 are satisfied. Then, as we have noted in (4), one can construct 
a process with law V{6i/2, . . . , On/'i) by an appropriate time-change of the 
BESQ processes. We continue the notation in (4). Recall the sum processes 

C = Er=i^iand5 = Er=i^i- 

Then, the market weights at any time t is given by the relation 

^'^^'^ = = J ^^^^^^ ' ^^^^ = 

Our first step is to get rid of the time-change by studying the process at 
random stopping time (,a when the process S hits a level o. 

Consider the corresponding hitting time (Ta for the process C,. Then, it 
follows plainly from the time-change relationship that S{t) = C(A(t)) that 
A(?a) = Ca and that 

(19) K^a) = 7 (A(?a)) = 7 K) = -Z [aa] ■ 

On the right side above we have the process Z the first time it escapes from 
the set a©. Since each 6i is positive the BESQ process can only exit © 
through the oblique boundary & (all the other boundaries are reflecting). 
Our objective is to compute this exit density, which in turn gives the exit 
density of the market weights at ?a . 

Before we head to computations, we remark that it is enough to take a = 
1. This is because of the following scaling property of the BESQ processes. 
Let y be a BESQ^ process. Then, for any positive a, the scaled process 
{a-^Y{at), t > 0} is a BESQ*^ starting from x/a. In particular, by scaling 
each Zi, . . . , Zn by a, we get that the law of the vector a~^Z(aa) is the 
same as the vector Z{ai) when the initial vector of values has been divided 
by a. 

The other consideration is whether 0"^ is flnite or not. The sum is a 
BESQ process of dimension Oq. This process is transient if and only if > 2. 
Thus, under the assumptions in Proposition 2, the finiteness of ai holds with 
probability one. 

2.1. Green kernel and the exit density of BESQ processes. Our main in- 
strument is to define a functional transformation analogous to the classical 
Kelvin transform. The intuition comes from the fact that when the dimen- 
sions of BESQ processes are positive integers, they have the same law as 
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that of the Euchdean norm-square of multidimensional Brownian motion. 
Thus the exit density from the unit simplex from the BESQ processes can, in 
principle, be derived from the Poisson kernel expansion for the exit density 
of the BM from the unit ball. One of the ways one can obtain the Poisson 
kernel formula is by employing the classical Kelvin transform techniques 
(see the book on Harmonic Function Theory [1, Chap 4]). We generalize 
that concept below. 

Consider the (generalized) Markovian generator of the process (Zi, Z2, ■ ■ ■ , Zn 
acting on (R""*"), the space of functions that are twice continuously differ- 
entiable in M""^ up to the boundary. It is the following differential operator: 

i=l i=l ' 

Any twice continuously differentiable function u that satisfies Cu = in an 
(open) domain D C will be called £-harmonic on D. 
Define the inversion map / : R"+\{0} — > M""*" by 

(21) I{z) 



It is easy to see that / is one-to-one and / o / is the identity map. Also, / 
inverts the interior of the punctured unit simplex 6\{0} to the interior of 
its complement in M""^\{0}. If is a domain in IR"'^\{0}, we will denote its 
image under the inversion map by I{D). 

Let D be a domain in M"'+\{0} and let u be a real valued function on D. 
One can define a function K\v^ : I{D) — ?• M as 

Y^^zA uiliz)), 00 > 2. 

Thus, K is a function that transforms a function on Z) to a corresponding one 
on I{D). We prove that it takes i3-harmonic functions on D to £-harmonic 
functions on I{D). We have the following proposition. 

Proposition 5. For any function u on D, define 

Then K[u] is a function on I{D), and we have 

/:K[u]{z) = K (z), for all z G 1{D). 
Thus, if u is C-harmonic, then so is K[u]. 
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To do the proof, we will need the following lemma. 

Lemma 6. Let p be a polynomial in n-variahles that is homogenous of 
degree m. Then, on any domain D C R"+\{0}, we have 

Proof. First note that for any two functions f,g, we have 

n 

(23) C{fg) = fC{g) + gC{f) + Zidifdig. 

i=l 

Now, for any power r we have 

, „ / \ / \ r-2 

n \ ' n I n \ n I n 



1 



r— 1 / \ — 1 

^E^'IE^^ I +2r(r-l) (f^z, 

\ r-1 



.4 = 1 



Thus, using the product formula (28), we get 



r-l 



(^E ^« j ^'(^) j = (^E j ^(^') + ^ (^0 + 2r - 2) [Y,^ z,j p{z) 

r-l 



n n 



+ 4^ E E 

j £(p) + r (00 + 2r - 2) ( p(z) + 4rm ( p(z). 

The final equality follows from the fact that for all homogenous polynomials 
of degree m, we should have (z, Vp) = mp. The easiest way to see this is to 
note that p{az) = a"^p{z) for all a > 0, take derivative with respect to a, 
and finally put a = 1. 
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Thus, we get 



r-l 



^[[^^ ^^^^ J " (^o + 2r-2 + 4m) j 

Choosing r such that + 2r — 2 + 4m = proves the Lemma. □ 



Proof of Proposition 5. We first prove this proposition when u is 
a polynomial p, homogeneous of degree m. By utilizing the property of 
homogeneity, we can write 



(24) 



CK[p] = C 



C 



^ n 

E 



i-eo/2 / 
P I 

l-6»o/2-2m 



p{z) 



l-9o/2-2m 



£(p). 



The final equality is due to Lemma 6. 

Now, note that, since p is homogenous of degree m, we have C{p) is 
homogenous of degree m — 1. Thus 



(25) 



K 



E 

\i=i 

' n 

E 



l-6»()/2-2 



l-6»(j/2-2-2(m-l) 



l-6»o/2-2m 



£(p). 



Combining equalities (24) and (25) we get 



liK{p\ = K 



which proves the proposition for the special case of homogenous polynomials. 

The general result now follows for all polynomials (obtained by taking lin- 
ear combinations of the homogeneous ones), and finally for all functions 
by taking suitable limits of polynomial sequences. □ 
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The explicit description of a Kelvin transform allows us to compute the 
Green function for the independent BESQ processes inside the unit simplex. 
As before, consider Z = be a vector of independent BESQ 

processes of dimensions . . . , 9^) respectively, satisfying the assumptions 
of Proposition 2. In that case, the process Z is transient (the sum Q being a 
BESQ process that is transient). 

Let (x, y) denote the transition density of BESQ^. Define the potential 
kernel of Z: 



Uy{x) = u{x,y) 



We compute this kernel below. 



Y[pt' {xi,yi)dt, x,ye 



i=l 



Proposition 7. The potential kernel u{x,y), when Oq = > 2, 

is given by the following formula. 



u{x, y) = n E r (^o/2 - 1 + 2m) 



s- 



-2m 



(26) 



i=l 



m=0 



ml 



m 



k: ki+...+k„=m ^ ^ i=\ ^ ^' 



n 



ixiyi) 



where S = J^l^^ixi + y^). 



Proof. The transition density of a BESQ process is explicitly described 
in [25, Appendix 7, p. 549] to be t~^f{y/t,6,x/t), where f{-,k,X) is the 
density of a noncentral chi-square distribution with /c-degrees of freedom 
and a noncentrality parameter value A. In particular, it can be written as a 
Poisson mixture of central chi-square (or, gamma) densities. Thus, we have 
the following expansion 



Pt 



{x,y) = t 



k=0 



where gr is the density of Gamma(r/2, 1/2). Taking products over 0j's we 
get 



(27) np*'(^-?/^)=*""n 



1=1 



i=l 



E 

k^=0 



-Xi/2t {Xi/2t) * 
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Recall the special notations from Subsection 1.1 to keep track of the 
various product terms. 

Since every term in (27) is nonnegative, we can expand the product as a 
series and get 

(28) ^ 

1=1 kl,...,kn *=1 

- E -"'^4p'>-"n™;7^^-"'''-''(7)"'"'"-"'=' 



L r(0i/2 + k,) 

k\ 



Here, j3{k) denotes the constant given by l/j3{k) = W^=i^{0i/2 + ki). 

To simplify this formula (28), it will be convenient to define S = S^+Sy 
^i{xi + yi). Thus, by regrouping terms, we get 



(29) n^'*'(^"y^) = y'^E«"'^'*(2*)"'"^'^"" E /5(^^)^ 

i=l m=0 ki+...+k„=m 

For notational convenience let us define 

k k 

(30) C(m)= E ^(^)^ 

ki+...+k„=m 

while we integrate out t from the expression in (29). 
Thus, we get 



poo ^ pc 

'■{x,y)= / T\pt'{xi,yi)dt = y''y2^(^) / 

JO ■ 1 n JO 

t=l rn=0 



Evaluating the inner integral is easy. Changing the variable to w = 1 /2t, 
we get 



oo roo 



-^/2*(2t)-^o/2-2m^^ 



1 f°° 1 

2 Jo 2 



2?i;2 

6*0 /2 
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Note that the assumption that > 2 is being used to show that the integral 
above is finite when m = 0. This completes the derivation of the formula 



^ oo k k 

u{x, y) = E r (^o/2 - 1 + 2m) 5-2^+^-^0/2 ^ /^(fe)^. 

m=0 fci>0, fci+...+fc„=m 

The expression in (26) can be obtained from above by dividing and multi- 
plying by m!'s inside the infinite sum. □ 

Lemma 8. The potential kernel Uy{x) is C harmonic in the interior of 
W^~^ . Moreover, for a fixed value of y, it has a uniform decay in the order 
OC^iXi)^~^°^'^ , as Xj tends to infinity. 

Proof. The first claim follows from the fact that u is the potential kernel, 
although it can be verified through direct computation. 

For the second claim, let Sx = ^^'^ ~ TliiVi- -'-'^ what follows we 

assume that Sx is much larger compared to Sy. 

Since we fix y and the vector 0, it follows from (26) that there is a constant 
C (depending on y and 9) such that 

(31) 



00 —2m / \ " 

m—n u- I — L — ^ \ i— 1 



ki\ 

m=0 k: ki + ...+kn=m i=l 

m=0 ^ ^ k: ki+...+kn=m ^ ^ i=l 

Recall the identity 



k: ki+...+k„=m i=l 

Since all terms considered are positive we get 

k: ki+...+k„=m i=l 
( 

n 



kj n E [k, kj n f = ^^^y 

fcl+...+fc„=m ^ i=l I I fci+...+fcn=m. ^ ^ i=l ' 
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Substituting this bound in (31) we get 

when Sx is large enough. This proves the claim. □ 

For the next proposition recall the stopping time cJi which is the hitting 
time of level one for the sum process C. 

Proposition 9. Consider the functional transform defined in (22), and 
define the kernel 

(32) v{x,y) = Uy{x) - K[uy]{x) = u{x,y) - K[uy]{x), x,yeM"+. 

Then, for every smooth nonnegative function f which is compactly supported 
away from the origin and any x £ &, we get 

(33) Ex r f{Z{s))ds= [ f{y)v{x,y)dy. 

Jo Je 

In other words, v is the Green potential on the unit simplex G for the process 
Z. 

Explicitly, the kernel v{x, y) is equal to 

(34) ^ 

m=0 



A;i>0: ki+...+kn=m 



ki--- knj f\ T{e,/2 + ki) 



Thus V satisfies the symmetry property 



n n 



(35) vix,y)llx'^'^'-' = v{y,x)lly','^'-\ 

1=1 i=l 

Proof. To prove this proposition we first note that for any compactly 
supported (in R""^) smooth test function / we have 

U{f){x):= f{y)u{x,y)dy = Ej f{Z{s))ds. 

JR"+ Jo 
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Thus, by the Markov property, it fohows that Mi(t) = U{f){Z{t))+ f{Z{s))ds 
is a martingale (^(0) = x), and that (Lemma 8) Cuy{x) = for all x ^ (3 
when y £ &. 

Fix an y G (3\(3. We now use Proposition 5 for the domain D = {x £ 
M""^ : Xi > 0, X ^ &}. It then follows that K[uy]{x) is /^-harmonic for all x 
in I{D), which is the interior of &. 



Claim. We now claim that, if we define 



W{f)ix) :-- 



f{y)K[uy]{x)dy, 



then M2{t) = W{f){Z{t A fii A a^)) is a martingale for every e > when 
Z(0) G 6. 

To prove this claim, it suffices to show that K[uy]{Z{t A cxi A cJe)), is a 
martingale for every y in &. We apply Ito's rule to the function A'[mj^]. Since 
K[uy]{x) is £-harmonic in the interior of ©, the process K[uy]{Z{tAaiAai:)) 
is a local martingale with the decomposition 

" d 

dK[uy]{Z{tAaiAae)) = 2^ ^y Zi{t A ai A a^) — K[uy]{Z{t))d(3i{tAaiAae). 



i=l 



The square bracket of this local martingale can be easily computed as 
d{K[uy]{Z)) (tAcJiAcr,) = 4^Zj(tAcJiAo-e) <^ —K[uy]{Z{t A ai A a,)) \ dt. 



Let us now compute the partial derivative: 



— K{Uy\{x) = — 

dxi ' ' dxi 



i-eo/2 



Xi 



U 



(1- V2) 



-9o/2 



U 



1-00/2 „ 



Here denotes the jth partial of Uy. 

Now, it can be seen from its explicit series expansion (26) that Uy is 
bounded and has bounded partial derivatives in when '^^yi < 1- Thus, 
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from the expression above, the partial derivatives of K[uy]{x) are uniformly 
bounded when x £ & and '^^Xi > e > 0. Hence it follows that K[uy]{Z){t /\ 
ai A (Te) is a martingale. By integrating with respect to f{y)dy, we have 
shown that M2{t Aa^) is a martingale for every e > 0. 

Thus the process 

N{t) = U{f){Z{t A cTi A a,)) - W{f){Z{t A di A a,)) + / / {Z{s)) ds, 

Jo 

is a martingale too. 

We now apply the Optional Sampling Theorem to this martingale at the 
stopping time aiAa^ (notice that the martingale is bounded by Ci+C2(cri Acre) 
which has a finite expectation). There are two cases to consider: when ai < 
(7e, we have Z{t A ai A a^) G S. For any x £ &, one has Uy{x) = K[uy]{x) 
and hence 

U{f){Z{aiAa,))-W{f){Z{aiAa,)) = when ai < a,. 
In the other case when < (Ji, one gets C(<7i A a^) = e. Thus 
(36) W{f){Z{ar A a,)) = e'~'"/^U{f){e~^ Z{ai A a,)). 

Thus, we get 

Ex r'''^ f{Z{s))ds= [ f{y)v{x,y)dy 
Jo Je 

+ Ex [-U{f){Z{a,)) + W{f){Z{a,)) \ a, < ai] P {a, < ui) . 

Since Q has dimension more than two, it almost surely does not hit the 
origin. Thus, it is clear that as e tends to zero, the left side of the above 
equation converges to Ex J^^ f {Z{s)) ds. We now show that the right side 
converges to Jq f{y)v{x,y)dy. 

Using the scale functions for it is easy to see that 

Now, as e tends to zero C/(/)(Z(cre)) remains bounded. The easiest way to 
see this is to note that / has a compact support away from the origin, and 
Z{ae) is away from all points in the support for a sufficiently small e. Hence 

lim Ex [-U{f){Z{a,)) \ a, < ai] P {a, < ai) = 0. 
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On the other hand, from (36) and (37) we get 
hm^, [W{f){Z{a,)) I a, < ai] P {a, < ai) = S^^/^ limE, [U{f){e~^ Z{a,)) \ a, < a^] 

The Hmit on the right is zero since €~'^Z(a^) tends to infinity and, thus, the 
function U (/) apphed to it uniformly goes to zero by Lemma 8. 
This completes the proof of the equality in (33). 

Now, we compute the kernel v{x,y) from the formula (26). Define a tem- 
porary notation 

_ r(go/2 - 1 + 2m) ^ f m \ A jxm)''' 

ml Ui---fcJiir(0,/2 + fc,)- 

Thus, from (26), we get 

i-eo/2 



v{x,y) = u{x,y) - \y2xi\ ^((^^0^'^) 



oo 



m=0 

u / -\ \ l-eo/2 oo / \ -2m 



m=0 

oo 

X; Rm{{S. + Sy)-'^+'-'°" - {S.Sy + l)-2-+l-^0/2| 



m=0 

This completes the derivation of the Green kernel. □ 

Finally, we derive the exit distribution from & for the process {Zi, Z2, ■ ■ ■ , Zn)- 
Note that the transition density can be guessed from the following version 
of Green's second identity for the generator of the BESQ processes. 

Lemma 10. Let C be the generator in (20). Let D C 6 be a compact 
domain with piecewise smooth boundary. Let m = (mi, . . . , m„) be the vector 
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given by mi = XiUi^x), where n is the outward unit normal vector at a 
boundary point x. 



Let ijj{x) be the weight function a;(x) := x^ = YYj=i ^. 



,72-1 



, and let u, v be 



two functions on D which are twice continuously differentiable ( continuous 
up to the boundary) on D. 

Then, assuming the right side below is integrable, we have 



(38) 



D 



dD 



{uCv — vCu) u}{x)dx = 2 
Here a is the surface Lebesgue measure on dD. 

Proof. Recall the generator C for BESQ processes: 



Cv 



d 



& 



d 



' dxi 



' dx 



dxi 



dv 
dxi 



We now use the Divergence theorem. Let n denote the outward unit nor- 
mal vector on a compact domain D with a piecewise smooth boundary. 
Given a vector-field of continuously differentiable functions F = (/i, . . . , /„) 
on D, define div F = Y^'^=idifi. Then 



div F[x)dx 



D 



dD 



F ■ n{y)a{dy). 



Here dD is the boundary of D, and (y{dy) is the surface Lebesgue measure 
on dD. 

Let /i be the measure on 6 given by the density oj{x). Now, we use the 
Divergence theorem to derive the following multivariate integration by parts. 



%/2 



dv 



du dv 



+ 



dD 



\ HQ ■ I I X ■ 

D~{ 'dxi dxi jj-^ ^ 
F ■ n{y)a{dy). 



dxi 

.,/2-l 



n 

dx 



The last equality above is obtained by applying the Divergence theorem to 
the function F = (/i, . . . , /«) where 



dv 

fi = n 4 



,72-1 



2n 



dv 
dxi 



11-7 



6^/2-1 



Interchanging u and v above and taking a difference, we arrive at (38). □ 
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Proof of Proposition 2. We first use the previous lemma to produce 
a convincing heuristic derivation of the transition density. 

Step 1: Heuristics. It suffices to prove the formula (7) for an arbitrary z 
in the open unit simplex. Fix any e > 0, small enough such that B{z,e) is 
contained in the interior of S. 

Consider a smooth nonnegative function / on (3. Let be a function on 
M"+ which is nonnegative, smooth, and zero outside B{z,€). We will use ip 
as an approximation to the delta mass at z. 

Consider the two functions defined on &: h{x) = J v{x,y)ip{y)dy, and 
l{x) = Exf {Zfj^). In the interior of 6, we have Ch = —ip by virtue of v 
being the Green potential, and l{x) is /^-harmonic as a corollary of its defi- 
nition. Assuming that both these functions are also smooth, with derivatives 
extending continuously to the boundary, we can apply the extended Green's 
identity (38) for u = I, v = h, and D = & to get 

(39) - / liy)i^iyMy)dy = 2 / (l^ - h-^) io{x)a{dx). 
JB{z,e) Jae V dm J 

Let us now analyze the right side of the above equation. The surface d& 
is piecewise linear and consists of the subsets Si, S2, ■ ■ ■ , Sn, and Sq = <B, 
where the outward normal vector for Si is — for i = 1,2, ... ,n, and for 
6, the vector is n~^/^l, the normalized vector of all ones. Thus, integrating 
separately on each Si, and dropping the constant 2 for the time being, we 
get 

(40) 

f (l^ - h-—] uj{x)a{dx) = / i^W- - ^TT^ ) w(x)cr((ix) = 
Jge\om dm J ~^Js,\om dm J 

- 1 (^1^ - ^1") -(-)-(^-) ^i^it-^ 0£ - ^1") ^^^^ 

Over each Si, for i = 1,2, ... ,n, the ith-coordinate Xi is zero. Due to the 
fact that each 9i > and assuming that h, I, and their partial derivatives 
are well-behaved, the integral above must be zero. Over S, by definition we 
have h = and / = /. Thus, combining (39) and (40), we get 

(41) - / l{y)i;{y)u;{y)dy = 2n-^l^ \ Y^^A /(^)7^ ^(x)n{dx). 
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Now, we take a sequence of '(/''s functions approximating the delta function, 
such that both h and its partial derivatives converge to and its corre- 
sponding partials. Thus we infer 

(42) l{z)uj{z) = -2 I f{x)uj{x) Xi—v{x, z)dx. 

The factor of l/-v/n gets absorbed when we parametrize the surface & by 
M"""^ (and hence dx represents the induced measure from the Lebesgue 
measure on M"~^). 

Since l{z) = Ezf{Z(j^), this identifies the exit density as 

a; G 6, z G 6. 

The problem with the above argument is that a priori we do not know the 
regularity of the exit distribution at the boundary of the simplex. But once 
we have guessed the solution, we can easily check that it must be the correct 
one. 



(43) Lpz{x) = ip{z,x) 



^uj{z) 



d 



.1=1 



x,z) 



Step 2: Computation based on heuristics. Let us now compute explic- 
itly the expression (43). To simplify matters, let us introduce some tempo- 
rary notations: for m = 0, 1, 2, . . ., let 



\i=l i-- 

E 

fcifciH — kn='m 



-2m+l-6»o/2 



Xi 



yi = l / \i = l J ) 



-2m+l-eo/2 



/ei • • • K.) n r(0i/2 + ki) ' 



{ZiX:-^^ 



Thus, from (32), we get 



d_ 

dxi 



z- ^ r(6lo/2 - 1 2m) 
2 ^„ m\ 



m=0 



n —n ^ R An 

aXi aXi 



imsart-aap ver. 2010/08/03 file: vsm_final.tex date: August 5, 2010 



24 



Now, when z is in the open unit simplex and x G S, we get 

/ n \ -2m-0o/2 

-^B^ = (-2m + 1 - dQ/o\ \ ..^^\ 



\i=l / 

_ (_2m + 1 - 00/2) ( 5^ zA ( + 1 j 

/ n \ / n \ -2m-6»o/2 

= (1 - 0o/2 - 2?n) ( 1 - 2^ + M 

We do not need to compute partial derivatives of since Bm is zero on 
6. 

Thus, by combining the partial derivatives of we get that ip{z,x) is 
equal to 



(44) 



m=0 \ 1=1 / \ i=l / 

which leads to formula (7) by substituting Sz for ^ - zi, noting that T{9q/2 — 
1 + 2m) (2m + 0o/2 — 1) = T{6q/2 + 2m) and completing the terms in the 
Dir density. 



Step 3: A rigorous proof. Now we show rigorously that the above formula 
is the true exit density. To do this we merely need to check that the heuristic 
derivation shown in Step 1 goes through. 

First we claim that 99(2:, x), as given by (43), is in the kernel of C in the first 
coordinate. That is, Cz'^{z, x) = for all z in the open unit simplex when 
X G (3. To see this, we use the symmetry property of the Green potential 
(35). Thus 



^{z,x) 



(45) 



^w(x) 
^uj{z) 



d 



d v{z,x) 



li=l 

n 



^ ^ dxi uj{x) 



d v{x,z) 

^ ' dxi u:{z) 
1=1 



Since Czv{z, x) = for all z in the interior of <3, it immediately follows that 
Cz^{z,x) must also be zero. 
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Now consider any smooth test function / on S, and following Step 1, for 
any y / in S\S, define 

(46) l{y)= lf{x)ip{y,x)dx. 

Je 

Now, from the explicit formula of ip{-, x) in (44) and usual analysis it follows 
that dip /dm (and hence 81/ dm) is a well-defined power series up to the 
boundary of ©. The function h in (39) is a convolution with a smooth 
mollifier ip, and is obviously smooth. Thus (39) and (40) go through and we 
arrive at the following modification of (42): 
(47) 

l{z)uj{z) = —2 ll{x)uj{x)y^Xi— — v{x,z)dx= ll{x)uj{z)ip{z,x)dx. 
Je ~[ oxi Jq 

We now cancel the factor ijj{z) appearing on both sides of (47), and com- 
pare with the definition of l{z) in (46) to get 



(48) I f{x)ip{z,x)dx = ll{x)(p{z,x)dx. 
Je Je 

Now we vary z in (3 and use the following claim. 

Claim. If li : 6 — 7- M is a bounded continuous function such that 

(49) lu{x)ip{z,x)dx = 0, for all z G 6, 
Je 

then u is identically zero. 

Assuming the claim for now, we get l{x) = f{x) for all x € &. Thus, the 
function I is harmonic with the correct boundary condition. The rest of the 
proof is now immediate. Since our density (p is supported over a compact 
set, it is enough to evaluate expectations of monomials under the density. 
For any monomial p = JlILi with each 7i > 1, we consider the function 



H{z) = lp{x)ip{z,x)dx, z £ &. 
Je 



le 

We claim that the process H{Z){t A o"! A a^) is a martingale for any e > 0. 
We have already shown CH = inside 6 (since Ctp = and taking the 
derivative inside the integral sign which is allowed since everything is smooth 
and bounded). Thus, the claim follows by noting from the explicit expansion 
of that the first partial derivatives of H are bounded on any domain 
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away from the origin. By applying the Optional Sampling Theorem we get 
H{Z{0)) = EH{Z[ai Acre)). One can now take e to zero, arguing exactly as 
in the proof of Proposition 9, to claim H{Z{Q)) = EH{Z{ai)). But, by (48) 
and the above claim, we see that H{Z{ai)) = p{Z{ai)). Hence we get 

E:,p{Z{ai)) = lp{x)(p{z,x)dx. 
J& 

Since the above identity holds for all monomials p, this completes the proof. 



It remains to prove the above claim. Note that since 93 is a power series, 
one can take the integral inside the sum in expression (7). Let (p^ be the 
resulting power series in z which is identically zero under (49). Let us define 
a change of variable qi = {\ + Sz)~'^Zi in expression (7) and expand ^p'^ as 
below: 

00 n 

(50) /'(z) = (l-5,)(l + S.)-^«/2 j;t(m) P(k)n9f' = 0- 

m=0 k\+...+kn=m i=l 

Here t{m) is the coefficient of the mth internal summand as in (7) and 

p(k) =(^^ j~ '"(^)^i^(^' ^ + ^/'^)dx- 

Since the right side of (50) is identically zero for all z £ &, it follows 
that the inner power series (as a function of the variables qi, ■ ■ ■ ,qn) must 
be zero. 

Now, fix any collection of positive integers fci, . . . , By taking repeated 
partial derivatives • • • d^"ip^{q), and letting each q^ tend to zero, we ob- 
tain p(k) must be zero for all k = {ki, . . . , /c„). But, from the structure of the 
Dirichlet densities, this shows that all multivariate moments of the measure 
u{x)uj{x)dx on 6 must be zero. Since © is compact this identifies u as the 
zero function. □ 

The proof of Proposition 1 follows immediately from combining Propo- 
sition 2, equation (19), and the discussion following it. One simply needs 
to keep in mind that in keeping with the time-change relationship (4), to 
compute the distribution of market weights under the model V(5i, . . . , (5„) 
we need to compute the exit density for n independent BESQ of dimensions 
9i = 2di for each i = 1, 2, . . . , n. 
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Readers might wonder if there is any direct way of seeing that the density 
expression in (7) (and hence (5)) integrates to one. Although not elementary, 
the following argument is such a direct way and serves as a sanity check. 

Assume, for notational simplicity, that 9q/2 is a positive integer r. Let s 
denote Sz- Since the integral of each Dir density appearing on the right side 
of (7) is one, by an application of Fubini-Tonnelli, we get 

(51) 

m=0 ^ fci+...+fc„=m ^ 

OO 

= (i-«)E 

m=0 

m.=0 ^ ^ 

OO 

= (1 - s)q ^ Po (§2m+r-l = 1 - r) . 
m=0 

Here refers to the law of a downward biased random walk {S^} with 
probability p of going up at each step and starting from a at time zero. 

Let Nk denote the (almost surely finite) number X]m=i ^{^m — then 
it follows easily that 

OO ^ 

V Po (S2m+r-l = 1 - r) = ^oA^l-r = 1 + Pl-r^Vl-r = 1 + ^oiVo = • 

^ — ' i — TT 

m=0 

Here vr is the return probability of the walk to zero, starting at zero. Since 
TT is well known to be 2p, one can substitute this value in (51) and get 

Je l-2p 
since s = p/q. This completes the verification. 

3. A skew-product decomposition result. The BESQ family of 
measures is well-known to be an additive family. This can be utilized to 
embed multidimensional BESQ processes in a measure-valued BESQ pro- 
cess as done by Shiga- Watanabe [26] and Pitman- Yor [24]. We follow the 
statement and notations in [24, Theorem 4.1]. 



2m + r-l\ (1 + ^)-2„.-.^^ 
\ m J 
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Let C[0, oo) be the canonical space of continuous paths with the usual 
topology. There exists a C[0, oo) valued process {Y^,0> 0,x > O) such that 
has law BESQ^. Moreover, we have the additive decomposition 

Y^ = Y^ + Yo', x>0, 9>0, 

where [Y^, x > 0) and (Yq , 9 > 0) are independent processes with sta- 
tionary independent increments, each of them having trajectories which are 
increasing and right continuous with left limits in C[0,oo). In other words 
both Yq and Y^ are independent C[0, oo)-valued Levy processes. 

Now, fix any nonnegative ^o- Let F be any distribution function (increas- 
ing, right continuous with left limits) on [0, ^o]- Consider the C[0, oo)- valued 
process (F^, < d < 6q), where 

^d = Y^^,) = Y^^^^+Y,^, 0<d<eo. 

Let 7^([0,6'o]) be defined as in Section 1.2 (iii). Given a realization of 
{Lrf, < d < ^o}; one can construct a P([0, 0o])-valued process n{t). For a 
fixed value of t and a sub-interval (a, b] in [0, ^o]) it assigns a mass 

fJ'{t){a,b\ = — — . 

This defines a probability measure uniquely, which we denote by the follow- 
ing notation 

Kt){A) = -1- r l{s G A)Tds{t), for aU A G B{[O,0o]). 
^ eai't) Jo 

We have the following skew-product decomposition result. Recall the defi- 
nition of the Fleming- Viot processes from Section 1.2. 

Proposition 11. Let co be the hitting time of zero for the process Tg^. 
Then, there exists a FV[0,6q] process {i'{t), t > 0}, independent o/Fgy, such 
that 

ft 

fi{t) = V (4Cf) , where Q = „ , t < o-q. 

Jo r0o(s) 

The proof is essentially one step away from the simpler finite-dimensional 
version that follows. 

Proposition 12. Let Z = {Zi, . . . , Zn) be a vector of n independent 
BESQ processes each of dimension Oi,. . . ,6n- Let C, be the sum 'Y^l^i Zi, 
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which is BESQ of dimension Oq = 9i + . . . + On- Assume that C(0) > and 
let 



(52) 



cTo = mf{t>0: C(t)=0}. 



Then, there is an n-dimensional diffusion v, independent of C,, and having 
law J(6'i/2, . . . , On/'i), for which 



(53) 



Z{t) = Q{t)v{ACt), Ct 



^ ds 



t < O-Q. 



Remark. The condition t < do is clearly necessary to guarantee that the time 
change Ct does not blow up. For n = 2 this result was noted by Warren and 
Yor in [29]. Also see the thesis by Goia [16]. 

Proof of Proposition 12. By our assumption each Zj satisfies the fol- 
lowing SDE: 



dZ.it) = e,dt + 2y^Z~{tjdl3,{t), i = l,2,...,n. 
Let Ri = ZijC,. Then the SDE for R^ for t < can be found by Ito's rule: 

dR,{t) = C^dZ,{t) + Z,{t)dQ-\t) + d{Zi,C^) 



= [Oidt + 2yz^dft(t)J + z,it) [-C^dCit) + C"^d(C> it)] - 4Z,C^ 
= [0iC^ - eoZiC^ + AZiC^ - AZiC^] dt 

n 

+ 2C^ ^/Zdf)dp,{t) - 2C^Zi{t) V^jdPj 

= C"' [0i - OoRi] dt + 2C'^/Zdt) [1 - C^Ziit)] dPi{t) - 2C^Z, 

n 

= [9^ - 9oR,] dt + C^/^2^/r,Y. (li^ = - \/^) 
Define the sequence of local martingales 

1/2 " 

(54) dMi{t) = J—— J2 = j} - /R^) d/3,{t), 
so that, 

(55) dRi{t) = [9, - 9oRi] dt + 2^ Ri{l - Ri)dMi{t). 
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However, since 

(1 - Rif + Ri^J = (1 - + (1 - ^0 R^ = {'^- Rt), 

we have guaranteed that (Mj) (t) = Ct- 

Let Tu be the inverse of the increasing function 4Ct, i.e., t„ = inf {t : Ct > u/A]. 
Let u = (z^i, 1^2, ■ ■ ■ 5 i^n) be the process obtained by time-changing R by r. 
In other words = Ri{Tu). Applying this time-change to the SDE for Ri 
in (55), we get 



(56) 



1 



dui{t) = - [9i - Ooui] dt + ^Jvi{l-Vi)Wi{t), 



where Wi is the Dambis-Dubins-Schwarz (DDS) Brownian motion (see [25, 
p. 181]) associated with Mi. This is the SDE for J{6i/2, . . . ,en/2) (see 
Subsection 1.2) once we prove that the diffusion matrix is given by a. To 
compute it, note that 



{ui,Uj){ACt) = {R^,RJ){t) 



Ri{t)Rj{t) = j] - ,/R~R, 



Now changing time by r, we immediately get {vi^Vj) = cr{i,j) as desired. 

All that remains now is to show is that the process above is independent 
of Q. The SDE for Q involves another martingale: 



(57) dC{t) = dodt + 2^Y^ ^/R'jdPjit) = 6odt + 2y^dp*{t). 

Here /3* is the local martingale J ^ ^jRjdfij , which is a standard Brownian 
motion by Levy's theorem ([25, p. 150]). Note that 



d(/3*,M,)(t) 



1 



Rill — Ri 



Thus, by Knight's theorem [25, page 183], the DDS Brownian motions of 
(Ml, . . . , Mn) and f3* are independent. This shows independence of {Wi, . . . , Wn 
and /3* . It is known ([25, page 439]) that (" is a strong solution of the SDE 
(57). Thus, from the independence proved above, it follows that C is inde- 
pendent of the vector {Wi, . . . , Wn), and hence v in (56). This completes the 
proof. □ 
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Proof of Proposition 3. It will be useful for us now to analyze the 
time-change Cj in Proposition 12. Let us define S{u) = C (t^), where r, used 
in the proof above, is the inverse of 4Cj. Since the derivative with respect 
to t of Ct is l/Cit), it follows that 

In other words, = S{t)dt. Thus, if we define Xi{u) = Zi{Tu), for 
i = 1, 2, . . . , n, it follows that 

1 r 

Xi{u) = Zi (tu) , Tu = - S{t)dt 

4 Jo 

which is exactly the solution of V{9i/2, . . . , On/2) as described in the Intro- 
duction. 

The first part of Proposition 3 is now established. The rest follows from 
known invariant distributions of the Wright-Fisher diffusions. See, for ex- 
ample, [11]. □ 

Proof of Proposition 11. Consider any finite sequence of Lebesgue 
measurable sets {Ai,A2, . . . ,An}. By our construction of the Levy process 
of BESQ processes, it follows that 

(58) Zi{t)= / l{se A,)rds{t), i = l,2,...,n, 

Jo 



are a sequence of independent BESQ processes of dimensions {6i, . . . ,6^ 
where 6i is the Lebesgue measure of Ai. Note that the sum ^ = Tg^ is a 
BESQ process of dimension 6q = 9i + . . . + On- 

By Proposition 12, there is a Wright-Fisher diffusion process i^(t) such 
that the time change relationship (53) holds for all t less than ctq. As before, 
let Tu be the inverse of the increasing continuous function 4Ct- 

One can then define a measure on the ci-algebra generated by { Ai , . . . , An} 
by defining 

u{u){Ai) := Vi{u) = —4 — zZiiTu), i = 1,2,... ,n. 

Note that in the pathwise construction (58), the time-change is the same for 
all choices of n and sets {^1,^25 • • • ) An}- Thus, it follows that the measure 
i^(ii) is consistently defined over any refinement of the sets Ai, . . . , An- More- 
over i/(n) is countably additive since it is derived from the measure TdsiTu)- 
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Thus, by a standard argument invoking the Caratheodory extension theo- 
rem, a unique probabihty measure is estabhshed on the Borel sets in 
[O,0o]- It is now clear that the entire measure-valued process u > 0} 

satisfies all the defining properties of the Fleming- Viot model as described 
in Section 1.2. □ 

3.1. Weights in a subset of the market. So far our analysis has considered 
the entire vector of market weights. It is often not possible to deal with all the 
stocks in a single large market. Transactions are expensive and the different 
market indices often concentrate on a chosen sub-collection of stocks. 

Thus, it is of interest to study the following problem. Suppose, without 
loss of generality, we consider the first m out of the total of n stocks in the 
equity market and define the process of sub-market weights as the vector 

X (t) 

(59) jl = ifli,...,flm), A.(i) = ^, S{t)=Y,Mt). 

^ ' i=l 

Can one describe the behavior of these sub-market weights ? The answer is 
yes, and the logic relies on a self-recursive property of the VSM models. Our 
next propositions makes this clear. 

Proposition 13. Consider the submarket weight vector ft as defined 
above. Then there exists a Wright-Fisher diffusion v, independent of the 
sum process X^iMi? such that 




du 



the equality holding for all t until fn hits zero. 



l,...,m, 



Since, we have already shown that the market weights have the same law 
as the Wright-Fisher models, we prove the Proposition for the WF models. 
We take, without loss of generality, m = n — 1, the case of a general m being 
on similar lines. 

Lemma 14. Let J = (Ji, . . . , J„) be the multidimensional diffusion J {Si, . . 
Consider the process 



Y 



J2 J3 J-n 



,1-Jl 1-Jl 1-Jl 

up to the stopping time ti = inf > : Ji{t) = 1}. 
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Then there is a diffusion v which is J{62,S3, . . . ,5n), independent of Ji, 
such that 

>-«-K/'r4(J))- °^*<- 

Proof. Let Zi,Z2, . . . ,Zn be n independent BESQ processes with re- 
spective dimensions given by 25i,262, • • • , 25n- Then, as was shown in Propo- 
sition 12, the process 

(60) Ji{t) = where = inf In > : 4 / ^>t], 

C{Tt) I Jo C(S) J 

is distributed as J {Si, . . . , 5n)- 

The proof utihzes the independence of the BESQ processes to derive the 
stated result. We first claim that the time-change ds/C{s) grows to infinity 
almost surely. To see this, note that, C is a BESQ process of dimension 
d = Yli=i ^i- When d < 2, the C process is recurrent, and thence the time 
change J"q (is/C(s) grows to infinity in finite time. 

When d > 2, it is known (see, e.g., [14, page 43]) that 

1 r ds 1 

(61) lim 



u^oo log u Jq ({s) d — 2 

Thus, again the time change grows to infinity with time u. The case when 
d = 2 can be sandwiched between the two cases above by using stochastic 
comparison theorems for BESQ processes. Thus, the process Ji(t)'s in (60) 
has been constructed for all time < t < oo. 

Now let Ci = '^^=2 sum of all the BESQ processes except the 

first one. Exactly as before, there exists a (n — 1) dimensional diffusion 
u = (^2, • • • , i'n), independent of Ci, with law J{62, ■ ■ ■ , Sn) such that 

/■* ds 

Zi{t) = Ci{tMAht), ht= / — -, i = 2,...,n. 

Jo Ci(s) 

Thus, for any i > 1, we get 

3-^^ = ||^ = ^a4MTt)), where h{n) = h^,. 
Let us now analyze the time change h{Tt). We get 



dt ' " ' " ' Ci(n) 1- Ji(t) 
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The computation of r/ has been done in the proof of Proposition 3. 
Thus, we get the description: for i = 2, . . . 



t 



tor ail r < Ti . 



Note that v is independent of both and Zi. Thus, P is also independent 
of Ji. This completes the proof of this result. □ 

4. Transition density of the market weights. Finally, we combine 
all the results we have derived so far to obtain the transition density for the 
market weights of the VSM model. 

Our first step is to analyze the stopping time in Proposition 1. To do 
this we return to the SDE (3) in the definition of the VSM model. As done 
in [13], one can express this SDE as 

5.. 



dXi{t) = -^S{t)dt + ^Xi{t)S{t)dWi{t), i = 1, 2, . . . , n. 

Summing up over all the coordinates, we recover the SDE for the process S 
as 

dS{t) = -S{t)dt + ^/WlY^ ^/ Xi{t)dWi{t) = -S{t)dt + 5(t)d/3(t), 

4 = 1 

where /3 is the local martingale / S~^^'^{t) ^ \/l(i{t)dWi{t), which is a stan- 
dard Brownian motion by Levy's theorem ([25, p. 150]). 

Thus, 5 is a Geometric Brownian motion, and can be alternatively ex- 
pressed as 

S{t) = 5(0) exp {{d - l)t/2 + /3(t)) . 

Thus, when S{0) = s, we get 

?a = inf {t > : P{t) + {d- l)t/2 > log(a/s)} . 

The density of the ?a is well-known and can be found in the Handbook 
by Borodin & Salminen [3]. To reduce clutter of notation, let us temporarily 
define 

j:={d—l)/2, and p = log{a/s). 

Note that p is assumed to be positive. 
Then the density of ?a given by 

- " exp(-<^)<i. 



V2^ V 2t 
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Recall now from Proposition 3 that the process S is independent of the 
market weights //. Thus, /x and are independent as well. Suppose we denote 
the transition density function of fi at time t by p{t,^,y), where the initial 
position ^ and the terminal position y are both elements of 6. Then, by the 
independence, it follows that 



(63) 



Jo 



where f is the exit density computed in Proposition 1. 

Note that, the above integral transform (63) has a single parameter p if we 
fix a = 1. The (63) is an integral transform. If one can invert this transform, 
we can recover p{t, y) from ip. As we show below, this integral transform 
is nothing but a Laplace transform in disguise. 

To see this, note that for any function h{t) (and keeping a = 1), we get 



1^°^ h{t)t-^/^e-^''/^exp (^-^^ dt. 



'2tt Jo 

If we now change the variable from t to u = we get 

(64) T(p) = / h{t)P (?i ^dt)=Pj= g{u)e-P ^'^du, 

Jo y Jo 

where the function g is defined by 

g{u) = /l(l/n)7x3/2e-7V(2n)^-2 ^ /,(l/^)^-l/2e-7V(2«). 

Thus, we get that 

(65) V2^Tip)p-'e~''^ = Aig) (pV2) , 

where A{g){-) represents the Laplace transform of g. In other words, the 
function g (hence h) can be recovered by inverting the Laplace transform. 

We are now going to apply the preceding analysis to the function h{t) = 
p{t, ^, y) for fixed values of ^ and y in the set S. In that case, from (63) and 
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the formula (5) (taking a = 1, x = e ^^), we get that 
Tip) = (1 - e-P) T ^^^"^ + (1 + e-")-'""-' 

m=0 ^ ' 



+...+kn=m 2=1 



X 

fe>0: fci+...+fc„=m 



Simphfying a bit, we get 

m!r(77T, + d) 



m=0 



+...+A;„=m j=l 



A;>0: A;i+...+A;„=m 

In particular, this nice series representation allows us to take the inverse 
Laplace transform inside the infinite sum and obtain the final formula. 

^ V{2m + d) 
^-^ mil [m + d) 



m=0 

X 



k: ki+...+kn=in i=l 

Here Qm{u) is defined by the Laplace transform formula 

A ie^) (pV2) = ^/2^p"^e-(™+^)'' (l - e~P) (l + e"'')-'"^-'^ , m = 0, 1, , 

Changing the variable back to t = we obtain 

hit) = Hl/u) = nV2e-^/(2«),(,) = f; ^g"^+/]. nV^e^V(2.),^(,) 

m=0 ^ ' 

\ ^ I m \ ~i r 

X 



fc: A:i+...+fc„=m i=l 

EJi^'-w E rite)'' D.fe;...). 

m=0 ^ ' k: ki+...+kn=m ^ ' i=l 



Here, the coefficients 6m 's are given by (see (64)) 



^"6„(t)t-3/2e-^'*/2exp dt = K{Qm){p'm- 



This establishes Proposition 4. 
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